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We study the statistics of the first passage of a random walker to ab- 
sorbing subsets of the boundary of compact domains in different spatial 
dimensions. We describe a novel diagnostic method to quantify the 
trajectory-to-trajectory fluctuations of the first passage, based on the 
distribution of the so-called uniformity index oj, measuring the similar- 
ity of the first passage times of two independent walkers starting at the 
same location. We show that the characteristic shape of P(oj) exhibits a 
transition from unimodal to bimodal, depending on the starting point of 
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the trajectories. From the study of different geometries in one, two and 
three dimensions, we conclude that this transition is a generic property 
of first passage phenomena in bounded domains. Our results show that, 
in general, the Mean First Passage Time (MFPT) is a meaningful char- 
acteristic measure of the first passage behaviour only when the Brownian 
walkers start sufficiently far from the absorbing boundary. Strikingly, in 
the opposite case, the first passage statistics exhibit large trajectory-to- 
trajectory fluctuations and the MFPT is not representative of the actual 
behaviour. 



1. Introduction 

The concept of first passage underlies diverse stochastic phenomena for 
which the crucial aspect is an event when some random variable of inter- 
est reaches a preset value for the first time. A few stray examples across 
disciplines include chemical reactionsp^ the firing of a neuron}^ random 
search of a mobile or an immobile target jSHUl diffusional disease spread- 
ing,^ DNA bubble breathingj 30 * 31 ! dynamics of molecular motorsJ^HaH the 
triggering of a stock option,^ etc. One distinguishes also between contin- 
uously varying variables for which the first passage across a given preset 
value coincides with the first arrival to exactly this value, and discontinuous 
processes for which this value can be overshoot, as it happens, e.g., for Levy 
flights characterised by long-tailed jump length distributions with diverg- 
ing variance! 36 * 37 ^ A variety of first passage time phenomena and different 
related results have been presented in Refs. I38|39l This volume presents an 
even broader exposition of relevant areas and describes the state of the art 
in physical and mathematical comprehension of such phenomena. 

In this chapter we will be concerned with first passage phenomena for 
particles executing Brownian motion. The distribution of first passage times 
for such processes occurring in unbounded domains is typically broad, such 
that not even the mean first passage time exists!^! In particular, in one- 
dimensional, semi-infinite domains the first passage time distribution of 
a Markovian process is universally dominated by the t -3 / 2 scaling nailed 
down by the Sparre Andersen theoremP3 A similar divergence of the mean 
first passage time occurs in stochastic processes characterised by scale-free 
distributions of waiting times ESES This signifies that sample-to-sample or 
trajectory-to-trajectory fluctuations are of a crucial importance for such 
processes so that the mean values — the mean first passage times, even if 
they exist — have little physical meaning. 

However, in many practically important situations first passage pro- 
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cesses involve Brownian particles moving in bounded domains (see, e.g., 
Refs. I42TI47]) . In this case the random variable of interest, e.g., the first 
passage time r to a boundary, a target chemical group, a binding site on 
the surface of the domain or elsewhere within the domain, etc., has a dis- 
tribution ^(r) which possesses moments of arbitrary, negative and positive 
order. Such distributions are usually considered narrow, as opposed to 
broad distributions, which do not possess all momentsP^D 

Nonetheless, even for such distributions first passage variables do fluctu- 
ate, of course, from sample to sample or from trajectory to trajectory, and 
an important question is how to quantify such fluctuations. One usually 
resorts to standard measures of the statistical analysis, such as, e.g., stan- 
dard deviations, the skewness or the kurtosis of the distribution. But these 
quantities are often not very instructive since they just produce numbers 
and it is not always evident how to interpret or compare them. Clearly, 
quantifying the trajectory-to-trajectory fluctuations in a robust and mean- 
ingful way is a challenging problem with many important applications. 

In this chapter we focus on this non-trivial problem and discuss a re- 
cently proposed procedure which presents a lucid illustration of the effect of 
the trajectory-to-trajectory fluctuations. It allows to quantify their impact 
and, moreover, shows that in some cases the distributions considered as 
narrow behave effectively as broad distributions 122123] ^ e em pi y a novel 
diagnostic method based on the concept of simultaneity of first passage 
events. Instead of the original first passage problem of quantifying the 
statistical outcome for a single Brownian walker, in this procedure one si- 
multaneously launches two identical, independent Brownian particles at the 
same position ro , namely two different realisations of a single Brownian Mo- 
tion (BM) starting at ro . The corresponding outcomes are the first passage 
times Ti and ti. One defines then the random variable 

n + r 2 

such that w ranges in the interval [0, 1]. The uniformity index ui measures 
the likelihood that both walkers arrive to the target simultaneously: when 
(J is close to 1/2, the process is uniform and the particles behave as if 
they were almost performing a Prussian Gleichschritt. In contrast, values 
of uj close to or 1 mean highly non-uniform behaviour, strongly affect- 
ing the trajectory-to-trajectory fluctuations. We note parenthetically that 
similar random variables have been used in the analysis of random probabil- 
ities induced by normalisation of self-similar Levy processes,^ of the fractal 
characterisation of Paretian Poisson processes,^ and of the so-called Match- 
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making paradox ] 50 ! 51 ! and more recently, to quantify sample-to-sample fluc- 
tuations in mathematical finances) 52 ! 53 ! chaotic systems the analysis of 
distributions of the diffusion coefficient of proteins diffusing along DNAs 55 
and optimal estimators of the diffusion coefficient of a single Brownian tra- 
jectory.^ 

To illustrate the concept of simultaneity, consider a generic, generalized 
inverse Gaussian form of ^/(r) (see e.g., the discussion in Rcf. 22 23] and 
references therein) 

^-exp^-l^exp^l), (2) 

where a and b are some constants. Note that the first passage time distri- 
bution in Eq. ([2]) is exact in the particular case of Brownian motion on a 
semi-infinite line in presence of a bias pointing towards the target site, or, 
equivalently, for the celebrated integrate-and-fire model of neuron firing by 
Gerstein and Mandelbrot.^ In general, the detailed form of vE'(t) for dif- 
fusion in bounded domains is obviously much more complex than given by 
Eq. ([2]), depending on the shape of the domain under consideration and the 
exact boundary value problem and is typically given in terms of an infinite 
series. 

Nonetheless, on a qualitative level Eq. ([2]) provides a clear picture of the 
actual behaviour of the first passage time distribution in bounded domains. 
Namely, ^(r) consists of three different parts: a singular decay for small 
values of r, which mirrors the fact that the first passage to some point start- 
ing from a distant position cannot occur instantaneously. This is followed 
at intermediate times by a generic power-law decay with exponent \x — the 
so-called persistence exponent) 57 ! 58 ! depending on the exact type of random 
motion and spatial dimension. Finally, an exponential decay at long r cuts 
off the power-law. A crucial aspect is that the exponential cutoffs at both 
short and long r ensure that in bounded domains ^(t) possesses moments 
of arbitrary positive or negative order. The parameters a and b depend on 
the shape of the domain, its typical size, and the starting position of the 
particle within the domain. When the linear size of the domain (say, the 
radius R of a circular or a spherical domain) diverges {i.e., R — > oo), the 
parameter b also diverges such that the long-time asymptotic behaviour of 
the first passage time distribution is of power-law form without a cutoff. In 
this case, at least some, if not all, of the moments of ^(r) diverge. 

Now, given the distribution ^(t), the distribution of the uniformity 
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index u> in Eq. ((TJ) can be readily calculated to giv d 53 * 54 * 

As we explicitly show in the following section [5J for the generic distribution 
in Eq. ([2| one finds the following explicit formP^l 



where K2 fl (-) is the modified Bessel function of the second type. 

One notices first that the form of the distribution P(ui) in Eq. ((4]) is 
distinctly sensitive to the value of the persistence exponent fi, which char- 
acterises the scaling behaviour of the first passage time distribution \l/(r) 
at intermediate times. Thus, for fi > 1, P(uj) is always a unimodal, bell- 
shaped function with a maximum at uj — 1/2, which implies that in such 
a situation both walkers (or two distinct trajectories of one walker) will 
most likely arrive to the target simultaneously. For fj, = 1, P(w) is almost 
uniform, P{u>) ~ 1, apart from narrow regions at the corners ui = and 
ui = 1, for b/a S> 1. Curiously, for fi < 1, which corresponds to the most 
common case, there exists a critical value p c of the ratio p = b/a such that 
for p > p c the distribution P(uj) has a characteristic M-shaped form with 
two maxima close to and 1, while at ui = 1/2 one finds a local minimum. 
Such a transition from a unimodal, bell-shaped to bimodal, M-shaped form 
mirrors a significant manifestation of sample-to-sample fluctuations. 

In what follows we survey known results and further explore this intrigu- 
ing behaviour of the first passage time distribution and of the corresponding 
distribution of the uniformity index ui for diffusion in bounded domains, fo- 
cusing on the effects of the domain shape, dimensionality, location of the 
target, the type of the boundary conditions, and of the initial position of 
the walker. 

2. First Passage Times and the Uniformity Distribution 

Consider a BM inside a general d-dimensional domain S, whose boundary 
dS = dS a U dS r comprises reflecting, dS r , and absorbing, dS a , parts. At 
time t = 0, the BM initiates at ro € S and evolves within the domain 
until the trajectory hits dS a for the first time at some random instant r. 
Furthermore, let P(r, t|ro) denote the conditional probability distribution 
for finding the Brownian walker at position r at time t, provided the initial 
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condition was at r at t — 0. The distribution P(r, t\r ) is the solution of 
the diffusion equation 

d t P(v,t\r ) = D d 2 r P(v,t\r ) (5) 

on <S, subject to the initial condition as well as the boundary conditions 
at dS. Here D is the diffusion coefficient and d 2 is the d-dimensional 
Laplacian. The solution of this initial boundary value problem is, in the 
best case, cumbersome, and explicit solutions may be obtained for only few 
simple geometries, (see e.o.Ref. |59|). 

If a finite part of the boundary is absorbing, i.e., dS a is not empty, then 
the distribution P(r, t\ro) is no longer normalised. The survival probability 
S? VQ (t) that the walker has not reached dS a up to time t, is defined by 

^r (t) = I P{v,t\v )dv. (6) 

Js 

^ ro (t) is a monotonically decreasing function of time, eventually reaching 
zero value, limt^oo J^ ro (t) = 0. In terms of the survival probability, the 
distribution of first passage times to the absorbing boundary becomes 

<Mr)=-^g^, (7) 

and the MFPT associated with the distribution ^(r) is defined as the first 
moment 



(r)(r ) 



poo pOO 

/ r* ro (r)dr= / y ra (r)dr. (8) 
Jo Jo 

In most of the existing literature, the dependence of the MFPT on the 
starting position of the walker is either simply neglected, or it is assumed 
that the starting point is randomly distributed within the domain S. How- 
ever, as we proceed to show, the ro-dependence of the first passage time 
distribution is a crucial aspect which cannot be neglected. 

From the First Passage Time Distribution (FPTD) one readily obtains 
the distribution of the uniformity index w defined in Eq. (TTJ), for two inde- 
pendent BM, in terms of its moment generating function 

$(A)=/ P(u) cxp (-Aw) dw, (9) 



o 



with A > 0. Since T\ and Ti are independent, identically distributed random 
variables, expression ([9]) can be formally represented as 

$(A)= / / ^(r 1 )^(r 2 )exp | -A - ) dr x dT 2 . (10) 

J0 J0 \ Tl+T 2 J 



May 6, 2013 2:10 



World Scientific Review Volume - 9in x 6in 



FPT-book 



First passage phenomena in bounded domains 7 

Integrating over dr\ we change the integration variable, t\ — > uj, so that 
Eq. (|10j) is rewritten in the form 



$(A) 



1 ( i \ duj 
exp(-Au;)^— 



t 2 *(t 2 )* 



1 -uj 



T2 d,T2- 



From comparison with Eq. ©, one obtains P(uj) in Eq. ([3]). 



(11) 



2.1. Trajectory-to-trajectory fluctuations and the shape of 
P(u>) 

To better grasp the relation between the trajectory-to-trajectory fluctu- 
ations and the distribution of the uniformity index, consider P{uj) of 
Eq. (@| corresponding to the FPTD of Eq. J2|). First, one readily no- 
tices that P(uj) vanishes exponentially fast when uj — > or uj —> 1 so that 
P(uj = 0) = P(uj = 1) = 0. Second, P(uj) is symmetric under the re- 
placement uj — > 1 — lu. Since the distribution vE'(t) possesses moments of 
arbitrary order, our first guess would be that P(uj) is always a bell-shaped 
function with a maximum at u = 1/2. 

To study the shape of the distribution of the uniformity index, we ex- 
pand P(uj) in Taylor series around its symmetric point w = 1/2 up to 
second order in (lu 



1) 4 

2) ' 



P(UJ) 



KUC) 



1 + 4 1 - u - C 



g 2M -i(2C) 



1 



(12) 



where C = Inspecting the sign of the coefficient before the 

quadratic term, i.e., 



g = l-(i-C 



K 2f <(2C) 



(13) 



we notice that 



For fj, > 1, g is always negative for any value of b/a so that here 
the distribution P{uj) is a bell-shaped function with a maximum at 
uj = 1/2. 

For a = 1, g is negative and approaches from below when b/a — > 
oo. It means that P{lo) is generally a bell-shaped function with a 
maximum at uj — 1/2, but it becomes progressively flatter when b/a 
is increased, so that ultimately P(uj) ~ 1 apart from very narrow 
regions at the edges for b/a^> 1. 
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Fig. 1. The distribution P(ui) of Eq. Jl): Panel (a), for fj, = 0.25 and different values 
of ft. Panel (ft), for b = 200 and different values of fi. Parameter a is set equal to 1. 



• For < fj, < 1 there always exists a critical value y c {n) of the 
parameter y = b/a which is defined implicitly as the solution of 
Eq. (|T3| for g = 0. For b/a < J/ c (m)> the distribution P(u>) is 
unimodal with a maximum at uj = 1/2. For b/a = y c (fJ-), the dis- 
tribution is nearly uniform except for narrow regions in the vicinity 
of the edges. Finally, which is quite surprising in view of the fact 
that in this case ^(r) possesses all moments, for b/a > y c {fj) the 
distribution is bimodal with a characteristic M-shaped form, 

two maxima close to and 1 and uj = 1/2 being the least probable 
value. 

We depict in Fig. [T] three characteristic forms of P(ui) for fj, = 1/2 and three 
different values of b/a. 

Therefore, for fi > 1 two independent BMs will most probably reach 
the target simultaneously. For < /i < 1 two distinctly different situations 
are possible: if b/a is less than some well-defined critical value y c (jJ<), then 
most likely both BM will arrive to the target for the first time together. 
If, on the contrary, b/a exceeds this critical value, the event in which both 
BMs arrive to the location of the target simultaneously will be the least 
probable one, yielding very large fluctuations in the MFPT obtained from 
different trajectories. 

In the rest of this chapter we present several concrete examples in which 
the transition between unimodal and bimodal distributions of P(uj) is ob- 
served. 
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3. One-dimensional Brownian Motions 

In this section we study the trajectory-to-trajectory fluctuations of the first 
passage time, in terms of the distribution of the uniformity index u for two 
independent identical BM. 

3.1. BM in the semi-infinite interval with a bias 

We first consider a BM in the semi-infinite interval S = [0, oo) with dS a — 
{0} and dS r — {oo}, in the presence of a constant bias pointing towards 
the target, for which the FPTD in Eq. ^ with fi — 1/2 is exact. In this 
case one has a = Xq/4D, where xq is the starting point, D the diffusion 
coefficient, and b = AD/v 2 (v being the drift velocity). Hence, the Peclet 
number, defined as the ratio between the rate of advection and the rate of 
diffusion, is Pe = 2^/ajb = xo\v\/2D (see, e.g., Ref . 131?)) . Consequently, we 
can make a following statement: 

Consider two independent identical BMs on a semi-infinite interval, 
starting at the same point x , having the same diffusion coefficient D and 
experiencing the same bias F which points towards the origin so that the 
drift velocity of both BMs is v < 0. Then, an event in which both BMs 
arrive for the first time to the origin simultaneously is 

- the least probable if Pe < Pe c , 
— the most probable if Pe > Pe c , 

where Pe c is the solution of the transcendental equation 

^^W^r- ™' (14) 

Therefore, the MFPT to the target, which in this case equals :co/M, might 
be an appropriate measure of the search efficiency for sufficiently large 
Peclet numbers, but definitely is not in case of small Pe. In the latter case 
the trajectory-to-trajectory fluctuations are significant and the mean value 
is not representative of the actual behaviour. 

3.2. BM on a finite interval 



We consider now BMs in the finite interval S = [0,L] with dS a = {0} and 
dS r = {L}, that are initiated at x — xq at time t = 0. In this case the 
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Fig. 2. Two BMs in the interval [0,L]. Panel (a): Distribution P(uj) in Eq. JT6]| for 
different values of x /L. Panel (b): FPTD f XQ (r) in Eq. (15) for L = 100, D = 1/2 and 
different values of xq/L. The dotted curve corresponds to a decay ~ r -3 / 2 . 



FPTD can be obtained exactly as the infinite sum 

2nD y, ^ fx ^ ^ 7r 2 (n + 1/2) 2 L>; 

n=0 



^o(^)--^E^(yJ eX P( ^ )' (l5) 



where the coefficients A n (^S-) = (n+ |) sin f 7r ( n +V 2 ) z o j Consequently, 
after some calculations the normalised distribution P(w) has the following 
fornix 

_ 2 _d_ y, sin (7r(m + l/2)f ) cosh (jm + 1/2)^1^(1 - f )) 
" ~^ w i^o m + 1 / 2 cosh (vr(m + 1/2) 

(16) 

Panel (a) of Fig. [2] shows the distribution P(uj) in Eq. (|T6]) for different 
values of the ratio xq/L, evidencing the transition in the shape of P(oj) at 
the critical value w 0.605. 

Furthermore, the moments of arbitrary order of vE'(t) in Eq. (I15p can 
be straightforwardly calculated. One obtains 

where -EizfcG) are the Euler polynomials. Consequently, the mean and the 
variance are given by 

M = -4Mfr)- < 18 » 
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and 




respectively. One can readily check that both statistical measures are mono- 
tonically increasing functions of xq / L and do not show any sign of a particu- 
lar behaviour at xq/L rs 0.605. The same occurs for higher order cumulants 
of the distribution *(r)P21 

This puts in evidence the superior sensitivity of the statistics of the 
uniformity index uj. As we show in sections 14.21 and 15.21 . the unimodal or 
bimodal shape of P(oj) classifies the parameter space into regions inside 
which the MFPT is either meaningful or meaningless. 



3.3. Arcsine law mechanism 

We found that any two BMs arrive to the target at progressively distinct 
times the closer they are initially to its location, and should most probably 
arrive together when they are far from it. Notwithstanding this counter- 
intuitive result, such a behavior is the same as the one behind the famous 
arcsine law for the distribution of the fraction of time spent by a random 
walker on a positive half-axisP^ 

Once one of the BMs goes away from the target, it finds it more difficult 
to return than to keep on going away. 

How does this relate to the properties of the distribution of the first 
passage time? In terms of the FPTD of Eq. (|15[) in the previous section, xq 
(a 1/2 ) and I (6 1/2 ) define the effective size of the region in which the decay 
of the FPTD is governed by the intermediate power-law tail. The larger 
this region is the larger the fluctuations in the first passage time become.^ 
This is shown in panel (b) of Fig. [5J All three curves show an exponential 
behavior for small and large values of r, as in Eq. @. For large r all three 
curves merge which signifies that at such values of t the characteristic decay 
time is dependent only on I. The lower cut-off is clearly dependent only on 
the starting point xq. Evidently the power-law region characterized by a 
decay ~ t~ 3 ^ 2 , grows with decreasing xo/l, i.e., when the BMs start closer 
to the target. 



May 6, 2013 2:10 World Scientific Review Volume - 9in x 6in FPT-book 



12 T. G. Mattos, C. Mejia- Monasterio, R. Metzler, G. Oshanin & G. Schehr 

4. Two-dimensional Brownian Motions 

In this section we analyse the statistics of ui for BM in different two- 
dimensional bounded domains. 

4.1. BM in a disc with a reflecting boundary 

We first consider BM in a disc of radius L, centered at the origin S = {r : 
|r| < L}, with dS a = {r : |r| = 1} and dS r = {r : |r| = L}, where r G K 2 
and I < L. Therefore, the target consist in a concentric disc of smaller 
radius I. 

Suppose next that a BM starts at some point at distance tq from the 
origin and hits the target for the first time at time moment r. The FPTD 
^r ( T ) is explicitly given by^l 

ID °° 

*r (r) = -y Mhro,L)exp(-%DT) , (20) 

n=0 

where Z is the normalization, 



A n (l,r ,L) 



\lL^(\ n L)-P{U^\ n l)) z 

dU (X n x) 



U (x) = Y (XJ)J (x) - J (X n l)Y (x); U^(XJ) 



dx 



=1 

while X n are the roots of the function Y (X n l)Ji(X n L) — J (X n l)Yi(X n L) , 
arranged in an ascending order, and Y n (-) are Bessel functions of the second 
kind. Note that X n depends on L and I. 

In this case the distribution of the uniformity index uj can be written 
aPl 

w-ttt*^^.^*.), (21) 

m=0 v ' 

where &(X) is the characteristic function of the first passage time distribu- 
tion in Eq. (HH defined by^l 



$(A) = 7 ^ \ — TtH ) ^\ / V \ ■ (22) 

Fig. |3] shows the FPTD in Eq. and P(u) in Eq. ^2TJ) for fixed L, 
D and r, and several values of rg. Here we observe the same trend as in 
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Fig. 3. Two BMs in the two-dimensional disc of radius L = 100 with a reflecting 
boundary and concentric absorbing disc of radius I = 5, and D = 1/2. Panel (a): 
Distribution P(ui) in Eq. i ]2 1 [ i for different values of ro/L. The dotted curve corresponds 
to P av (oj). Panel (b): FPTD *,.q(t) in Eq. J20} for different values of r /L. The dotted 
curve corresponds to a decay ~ 1/rlog 2 t. 

the previous one-dimensional example, namely the FPTD becomes broader 
when the starting point of the BMs is closer to dS a , with an intermediate 
power-law behavior with a logarithmic correction ~ l/rlog 2 (r) (see the 
dashed curve in Fig.[3]-&). The critical value for the transition between bell- 
and M-shaped P(uj) is numerically found to be tq/L ~ 0.643. Moreover, if 
the starting point of two BM is averaged over the domain S, one readily 
finds, by averaging the result in Eq. (|2Tj) . that P av (w) = 1 (see the dashed 
curve in Fig. [3^6) . 

4.2. BM in a disc with aperture 

For other simple geometries but mixed absorbing/reflecting boundary con- 
ditions, it is not possible in general to obtain analytical closed expressions 
for the FPTD. However, the examples discussed previously suggest that the 
arcsine law mechanism in section 13.31 underlying the appearance of large 
trajectory-to-trajectory fluctuations of the first passage time, is generic in 
bounded domains of arbitrary shape. 

For instance, consider a circular domain of unit radius (in polar coordi- 
nates r = (r, 6)), S — {r : r < 1}, and the following boundary conditions: 
the segment with \9\ < 0/2 is absorbing while the remaining part of the 
outer circle is reflective. The aperture of the circular domain is thus an arc 
of length 0. 

Obtaining an expression for the FPTD involves the cumbersome prob- 
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Fig. 4. Phase-chart for the shape of the uniformity distribution P(lo) for a BM in the 
unit disc with reflective boundary (solid line) and an aperture of size corresponding 
to the absorbing boundary (dashed line), (a) = 7r/18, (b) = it/2, (c) = it, and 
(d) = 2-7T. Starting locations are coloured light blue if \ < — \+, dark blue if \ > X* 
and beige if |x < X*i where \* = !• 



lem of finding the solution of the diffusion equation with mixed Dirichlet- 
Neumann boundary conditions. Here we content ourselves with the numer- 
ical analysis of the problem. Let each BM commence at position (po,Oo) 
inside the unit circle and determine the set {r^} of first passage times to 
the location of the aperture. From these data one obtains P{oj)- To easily 
grasp the geometric dependence of the first passage time fluctuations, we 
considered the parameter \ that quantifies the shape of the distribution 
of a;p2l Performing a fit of the numerically obtained P(lu) to a quadratic 
polynomial of u> in the domain 0.05 < ui < 0.95, we define x as the coeffi- 
cient of the quadratic term, so that the sign of \ determines the shape of 
P(cli): x < corresponds to the unimodal, bell-shaped distribution, \ > 
signifies that the distribution is bimodal, M-shaped, and a zero value of 
X = means that P(w) is uniform. 

In Fig. |4] we show the phase-chart for the shape of P{oj) for different 
sizes of the aperture @- As in all previous examples, the closer the starting 
position of the BM is to the absorbing boundary the larger the fluctuations 



a Note that the case shown in Fig. [4] (d) reduces to a one-dimensional problemPH 



May 6, 2013 2:10 World Scientific Review Volume - 9in x 6in FPT-book 



First passage phenomena in bounded domains 15 

of the first passage time become. One thus expects that very near the ab- 
sorbing boundary the standard deviation of the first passage time becomes 
much larger than its mean. We note that this result is generic irrespectively 
of the aperture size. 

The same analysis was recently carried out for other different shapes, 
such as pie- wedges or triangles in Ref. H21 with the same conclusions. 

5. Three-dimensional Brownian Motions 

5.1. BM in a sphere with a reflecting boundary 

Consider a BM in a sphere of radius L (in spherical coordinates r = 
(p, 8, </>)), S — {r : p < L} with reflecting boundary dS r — {r : p — L}, and 
a concentric absorbing sphere of radius I < L, dS a — {r : p = I}- 

For this geometry, the distribution of the first passage time r of a BM 
starting at a distance po from the origin to the surface of the target is given 
explicitly bjJ^H 

D °° 

*po(t) = z22Ml,Po,L)exp(-\ 2 n DT) , (23) 

?i=0 

where the coefficients 

2u Q (X n p )u' (X n l) 
A n (l )P0 ,L)= g(JjjMt0 , 

with iio(x) = yo(X n l)jo(x) - jo(X n l)yo(x), u' (X n l) = duo(X n x)/dx\ x =i, and 

G(l,L,X n ) = Rj (X n L)u (X n L)(X n Ljo(X n l) - y (X n l)) - 
_ jo(X n l)yo(X n l) L - I 

X n l 2 ^n 

The set {A n } are the roots of yo(X n l)ji(X n L) — jo(X n l)yi(X n L), arranged in 
an ascending order, while j n (-) and y n (-) are the spherical Bcssel functions 
of the first and of the second kind, respectively. 

The corresponding distribution P{uj) for two independent BMs starting 
at a distance po from the origin ia^ 



oo 



n.m— 



A n (l, po,L)A m {l, p , L) 
(uXl + (1-^X1) 



2 



lyA m (lpo,L)^/ = l_u m2 , 



Z ^ XI 

m=0 
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Fig. 5. BM in a sphere of radius L = 100 with a reflecting boundary and concentric 
absorbing sphere of radius 1 = 5, and D = 1/2. Panel (a): Distribution P(to) in Eq. (124[ i 
for different values of po/L. The dotted curve corresponds to P av (o;). Panel (ft): FPTD 
>&p (t) in Eq. l|23|l for different values of po/L. The dotted curve corresponds to a decay 



where $(A) is the characteristic function of 4'(t) in Eq. ([2"3")) . defined b; 



sinh 



$(A) 



4 i? cosh \/%(L- po 



Po 



sinh 



.(L-0 



i? cosh 



(L-0 



(25) 



The distributions in Eqs. (|23|) and ((24)) are shown in Fig. [5] for fixed 
L, I and 13, and several values of po. An intermediate power-law ~ r~ 3 / 2 
is apparent for po/L = 0.055, persists for one decade for po/L = 0.3 and 
is entirely absent for po/L = 0.9, i.e., when the BM initiates close to the 
reflecting boundary (see Fig. [5]-&). Correspondingly we find that P(w) has 
a different modality depending on the value of the ratio po / L with a critical 
value of po/L ~ 0.68 (see Fig. [5]-a). For p /L > 0.68, the distribution has, 
in principle, a maximum at ui = 1/2 but this maximum is almost invisible 
so that visually the distribution looks more like a uniform one, as compared 
to the ID case in which the maximum is more apparent. When the starting 
point po is uniformly distributed in S, we again find P av (uj) = 1 (dashed 
curve in Fig. [5}a). 



5.2. Narrow Escape Problem 

As a last example we consider the problem of the so-called Narrow Escape 
Time (NET). This problem, generic in cellular biochemistry, consists of es- 
timating the time that a randomly moving particle spends inside a bounded 
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Fig. 6. Cut of the phase-chart for the shape of P(ui) for the NET problem with p a = 0, 
a = 1/4tt and D v = 0.000416 = D s = 0.000625. Starting locations are coloured light 
blue if x < ~X*> dark blue if \ > X* an d beige if \x\ < X*i where \* = 0.1- The red 
region on the sphere indicates the position of the narrow aperture. The phase-chart has 
rotation symmetry with respect to the position of the aperture. 

domain before it escapes through a narrow aperture on the boundary of the 
domain. The particle can be an ion, a ligand or a molecule diffusing inside 
a cell, a microvesicle, a compartment, an endosome, a caveola, a dendritic 
spine, etc., and a variety of processes in which the importance of the NET 
problem is striking have been discussed in the past.^^H 

In section 14.21 we discussed a two-dimensional version of this problem. 
In three dimensions, for a spherical domain of radius L, reflecting every- 
where except at a narrow circular aperture of radius a <C L in which the 
boundary is absorbing, the survival probability at large times is <5?r {t) ~ 
exp(— t/rsu), where the characteristic time decay tsu = irL 3 /'3Da corre- 
sponds to the MFPTP 

Here we discuss the statistics of the uniformity index for the NET prob- 
lem. Consider a sphere of unit radius S = {r : p < 1} (in spherical coor- 
dinates r = (p, 9, (/))), and the following boundary conditions: a solid angle 
| Q | < a is absorbing while the remaining part of the sphere is reflective. 
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Furthermore, consider a BM inside S with a diffusion coefficient Dy, and 
initially located at (jjq, #o, <A))- When the Brownian particle hits the sphere 
dS it becomes weakly absorbed with probability p a and starts to diffuse, 
with diffusion coefficient Ds, along the surface of the sphere. During a time 
interval At, the particle detaches from the sphere with probability p a At, 
and continue its motions until it hits the narrow aperture at time r. 

Considering couples of BMs one obtains P(w) as before. Fig. [5] shows 
the phase-chart for the shape of P(w), as defined in section l4~2l for p a = 0, 
i.e., when the BM never gets attached to the sphere. The results qualita- 
tively resemble those obtained in section 14.21 for the two-dimensional disc 
for narrow apertures (see Fig. @Ja). Note that the statistics of the first 
passage time inherit the polar symmetry of the geometry with respect to 
the position of the aperture. 

Our results indicate that for BMs starting near the location of the aper- 
ture, the MFPT has little significance to determine the NET. We have also 
considered finite values of the attachment probability p a (finite particle- 
surface affinity), obtaining qualitatively similar results than the ones shown 
in Fig. [6] for p a = 0. 

6. Conclusions 

In this Chapter we discussed a novel diagnostic method which allows to 
quantify straightforwardly and meaningfully the effects of the trajectory- 
to-trajectory fluctuations in first passage phenomena characterised by "nar- 
row" distributions 4 r (r) possessing moments of arbitrary order. This diag- 
nostic is based on the simultaneity concept of the first passage events. To 
this purpose, instead of the original random variable r characterising the 
first passage event, we introduced a variable lj = 7i/(ti +T2), where t\ and 
T2 are identical independent random variables with the same distribution 
\P(t). Physically, the uniformity index w shows us the likelihood of the 
event that two different trajectories arrive to the prescribed location for 
the first time simultaneously. 

We demonstrated that the characteristic shapes of the associated dis- 
tribution of the uniformity index to are very sensitive to the trajectory- 
to-trajectory fluctuation. In some cases, when such fluctuations have only 
little effect so that the MFPT can be considered as a meaningful property, 
the distribution P{uj) is unimodal with a maximum at ui = 1/2, such that 
T\ = T2 is the most probable event. In other cases, P(u>) is shown to have 
a bimodal form with two maxima close to lj — and uj = 1 and a local 
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minimum at u) = 1/2. This means that in such situations two different first 
passage events wili be most likely characterised by very different n and T2 , 
so that the trajectory-to-trajectory fluctuations will be very important and 
the MFPT cannot be considered as a meaningful characteristic measure 
of the process. In this regard, the underlying "narrow" distributions *(r) 
may show a behaviour specific to broad distributions which do not possess 
all moments. 

To illustrate this concept, we considered several examples of Brownian 
motion in bounded domains, in different spatial dimensions and with dif- 
ferent locations of the target. We demonstrated that, remarkably, within 
the same domain the very shape of the distribution P{lo) depends crucially 
on the location of the target and on the starting point of the Brownian 
motion — trajectories starting from one point may exhibit considerable 
fluctuations with respect to the first passage to the target, while trajecto- 
ries starting from some other point may arrive to the target location almost 
simultaneously. For some geometries, we evaluated the "phase charts" for 
the distribution P(co) showing which shape — unimodal or bimodal — this 
distribution will have for different starting point within each domain. 

We expect that a similar behaviour will be observed in bounded one- 
dimensional geometries for fractional BM with arbitrary Hurst index H or 
for a-stable Levy flights with < a < 1, and more generally, for bounded 
geometries with compact or non compact exploration. 
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